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INTRODUCTION 

Shock-wave unsteadiness has long been observed in rapidly compressed supersonic 
turbulent boundary-layer flows with significant separation. 1 The shock oscillations 
are characterized by large time scales, relative to those associated with the turbu- 
lence in the upstream, undisturbed boundary layer, and have been identified with 
low frequency, large amplitude (again, relative to the upstream boundary layer), wall- 
pressure fluctuations. 2,3 Researchers have been interested in these flows for some time 
because of direct aerodynamic and turbo-machinery applications. However, until re- 
cently, most experiments reported long-time-averaged data without addressing the un- 
steady aspects. This was often due to a lack of the necessary instrumentation. 

In the last few years, the advent of miniaturized high-frequency pressure transducers, 
and the continued development of the laser Doppler velocimeter (LDV), which non- 
intrusively measures velocity, have given new, timely, facility to tne study of shock- 
wave unsteadiness as it relates to boundary-layer behavior. This comes as demands 
for more detailed information on separated shock boundary-layer interactions have 
been generated by expanding efforts in computational fluid dynamics (CFD). To date, 
computations have been successful in predicting a limited class of such interactions, 4 ’ 5 
principally 2-D flows, and 3-D flows in which the shock waves were planar. The same 
combinations of numerics and turbulence models have as yet failed to predict a whole 
other group of separated flows, to be distussed here, with non-planar shocks. 6 ’ 7 This 
failure, coupled with knowledge that the turbulence models employed do not provide 
for the shock unsteadiness which is observed experimentally, has produced interest in 
the effect that unsteadiness has on the separating boundary layer. 

Certain basic questions have been raised. First, can the effects of the unsteadiness 
be removed from the mean experimental data by some conditional sampling tech- 
nique, so that “snapshots* of the flow, corresponding to specific shock positions (or 
phase angles), can be extracted and studied independently? Second, are the mea- 
sured turbulence quantities for flows with large separation zones dominated by some 
“pseudo-turbulence* which results from the overall bimodal nature (shock forward 
- shock back) of the interaction? Third, does it really matter? Can the unsampled 
time-averaged data be meaningfully compared to the computations, ignoring the fact 
that the shock wave is unsteady? If not, can a shock unsteadiness model be developed 
that would improve the predictive accuracy of current computational schemes? 


In recent years, a joint experimental and computational program aimed at address- 
ing some of these questions has been in progress at the NASA Ames Research Cen- 
ter. A Mach 2.85 shock-wave/turbulent boundary-layer flow was set up over a series 
of cylinder-flare bodies in the High Reynolds Number Channel I. In the first phase of 
this investigation, the transition from fully attached to fully separated flow was stud- 



ied using axisymmetric flares with increasing (12.5°-30°) compression angles. 8,9 In the 
second phase, the 30° flare was inclined relative to the cylinder axis, so that the effect 
on a separated flow of increasing 3-dimensionality could be observed. 6,7 For each of 
the separated cases investigated in phases 1 and 2, shock unsteadiness was discussed, 
but no attempt was made to remove its effect from the reported mean data. 

The current paper examines in some detail two 3-D separated cases. A simple condi- 
tional sampling technique is applied to the data to group them according to an asso- 
ciated shock position. Mean velocities and turbulent kinetic energies, computed from 
the conditionally sampled data, are compared to those from the unsorted data and to 
computed values. (The predictions are included to inform the reader as to the present 
state of the art. They do not represent any previously unpublished developments in 
CFD.) Finally, of the “basic questions” raised above, the first three are addressed. 
Namely, 1) can conditional sampling be used to provide “snapshots” of the flow; 2) are 
averaged turbulence quantities dominated by the bimodal nature of the interaction; 
and 3) is the shock unsteadiness really important to computational accuracy. 

DETAILS OF THE EXPERIMENT 

Complete details of the experimental setup and procedures have been previously 
reported. 6 A summary of the relevant aspects follows. 

The experiment was performed in a 25.4-by-38.1 cm rectangular test section of Ames 
Research Center’s High Reynolds Number Channel-I. This air-charged blowdown tun- 
nel employs interchangeable nozzles and test sections to generate a range of transonic 
and supersonic Mach numbers. 


CONE AXIS 



Figure 1 shows the basic flow configuration, and the associated coordinate system 
used in the data analysis. An axisymmetric, equilibrium turbulent boundary layer 
forms over the 5.08-cm-diam stainless steel cylinder which is aligned with the tun- 
nel axis. The cusped nose prevents formation of strong shocks that might otherwise 
reflect off tunnel walls back to the cylinder. The flow is abruptly compressed as it 
approaches the 30° half-angle flare, and the boundary layer then separates from the 
cylinder. For nonzero flare-inclination angles (a), the interaction is three dimensional. 


The principal region of study extended approximately five undisturbed boundary-layer 
thicknesses, do, upstream of the compression corner, and a similar distance down- 
stream along the flare. This paper discusses two 3-D configurations, a = 5° and 23°. 
The data presented are for the upper symmetry plane, z = (f> = 0, only. 

Nominal total and undisturbed flow conditions were: X< = 270° K, p t = 1.7 atm, 

Mqo = 2.85, Re = 16 x 10 6 /m, and do = 1.1 cm. 

Velocity measurements were made using a two-component laser Doppler velocimeter 
(LDV) operating in the forward scatter mode. The system employed a 4-W argon-ion 
laser, signal-processor counters, Bragg cells for frequency shifting, and half-micron la- 
tex light-scattering particles. Two-color operation of the LDV was used to distinguish 
the individual velocity components, while directional discrimination for each compo- 
nent was achieved with frequency shifting. Reference 6 gives details of the system 
components, and the procedures for alignment, calibration, and operation. A PDP 
11/34 minicomputer was used for data acquisition and initial analysis. To increase 
speed, the data were stored on hard disk during tunnel runs. They were later trans- 
ferred to magnetic tape for reduction, and then to a VAX 11/785 computer for all 
subsequent analyses. 

The cylinder portion of the model was instrumented with six high-speed pressure 
transducers (DC to 100 KHz) mounted just below the surface, and spaced 0.5 cm 
apart along the <b = 0° surface line. They were used during LDV data acquisition 
to record instantaneous wall pressures in the vicinity of the separation shock. 

DETAILS OF COMPUTATIONS 

The equations used to describe the mean flow were the time dependent, Reynolds- 
averaged Navier-Stokes equations for a 3-D fluid. For turbulence closure, the two 

equation, fc-e eddy-viscosity model 10 with wall-function boundary conditions 4 was 
used. The MacCormack explicit, second order, predictor corrector, finite-volume nu- 
merical procedure 11 was employed. 

The computational domain consisted of uniform mesh spacing in the streamwise (z) 
and azimuthal ( <j > ) directions. In the vertical (y) direction, a geometrically stretched 
spacing was used near the solid surfaces, followed by a uniform spacing. The total 
mesh size was 64 points in the streamwise direction, 33 in the vertical direction, and 
38 in the azimuthal. Typically, 16 mesh points were used to resolve the boundary 
layer. 

The upstream boundary conditions were prescribed by a combination of uniform free- 
stream conditions and the results of a boundary-layer computation matching the ex- 
perimental displacement thickness. At the downstream boundary, the gradients of the 
flow variables in the streamwise direction were set to zero. In the azimuthal direction, 
symmetry conditions were applied at <f> — 0° and <f> = 180°. No-slip and constant 
wall temperature conditions were imposed at the surface where turbulent kinetic en- 
ergy and dissipation were set to zero. Free-stream conditions were applied at the outer 
boundary. 

To reduce the time required to achieve a converged solution, wall-function bound- 
ary conditions were used at the cylinder and flare surfaces. The compressible two- 
dimensional wall functions derived in ref. 4 were extended to 3-D flows by replacing 
the horizontal velocity with the total velocity parallel to the wall. The flow direction 
was also assumed to remain constant between the surface and the first grid point away 
from the surface. The y + location of the first grid point varied from 30 to 120. 

CONDITIONAL SAMPLING 

To examine the effects of shock unsteadiness on the time-averaged data, such as veloc- 
ities and turbulent stresses, it was desired to construct different ensembles of instanta- 
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a) Shock forward. 



b) Shock back. 


Figure 2 Shadowgraph and related instantaneous surface pressures. 


neous measurements, each corresponding to a particular shock wave position (or phase 
angle). Then, by comparing the mean flow quantities derived from these ensembles 
with each other, and with those derived from the ensemble which embodied all shock 
positions, the response to changes in shock location could be observed. 






To accomplish this, it was necessary to conditionally sample, and sort, the instanta- 
neous LDV measurements according to the position (or region) the shock was in when 
they were acquired. This could be done during data acquisition itself by employing a 
trigger that, based on a pre-set condition, would have the LDV selectively take data. 
However, it was judged a more practical approach to conditionally sample during 
post-test analysis of the measurements that were made unconditionally, i.e., whenever 
a valid LDV event was registered. This required large numbers of data to be taken so 
that each sorted ensemble was large enough to yield statistically valid mean results. 
The question of how to link individual data with corresponding shock positions also 
had to be answered. 

Figure 2 illustrates the experimental procedure used to relate the shock movement to 
surface pressures. Simultaneous measurements were taken of the shock wave motion 
(with high speed — 9000 frames/sec — shadowgraph movies) and surface pressures at 6 
axial locations 0.5 cm apart (using high frequency miniaturized transducers). Com- 
puter generated timing marks were recorded on the film when surface pressures were 
sampled, and the result was that for each film frame there was an associated instanta- 
neous pressure distribution. The figure shows a high speed shadowgraph frame and its 
related pressure distribution for a “shock forward” and a “shock back” case. Differen- 
tial pressure transducers were used, and the values shown represent changes from the 
undisturbed wall pressure. 
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Figure 3 Correlation of shock position with surface pressure. 

Figure 2 also illustrates the definitions of X e , the streamwise distance from the cor- 
ner to where the shock intersects the boundary layer edge, and X p , to where the dif- 
ferential pressure crosses a value of 0.5. The variations of X 4 and X p with time, for 

a typical 60 msec interval, are shown in figure 3. The high degree of correlation be- 
tween these two quantities — they differ, for the most part, by a constant — indicates 
that shock location can be deduced, with reasonable accuracy and relative ease, from 
surface pressure measurements. 

All LDV data were measured and recorded simultaneously with the pressures ob- 
tained from the array of transducers. Numerous criteria, applied to the pressures, 
were considered for conditional sampling of the instantaneous velocities. However, for 
the present report, one of the simpler algorithms was chosen. 

For each LDV measurement location, the overall mean and standard deviation of the 
pressure (p and p' rmB ) were computed from the transducer which was positioned be- 
neath the average shock position, as determined visually from shadowgraph movies. 
Then, the instantaneous pressure levels (p(t)) were examined at that transducer, and, 
based on them, the corresponding velocity pairs (u(f), t/(f)) were selected for, or ex- 



eluded from, conditionally-sampled ensembles. This paper identifies a “shock forward” 
(SF) case which includes all velocity pairs for which p(t) > p + p' rm at the surface 
location beneath the mean shock position. It also presents a “shock back” (SB) case 
for which velocities correspond to p(t) < p - .5 p' rma . Finally, the flow quantities 
derived from the full data ensemble with no conditions applied, constitute the “all 
data” (AD) case. The coefficient of 0.5 in the SB criterion was necessary to retain a 
sufficient number of data to yield reliable statistics. The nature of the distribution of 
pressure levels is such that the number of measurements corresponding to SB drops 
rapidly as that coefficient is lowered below —0.5. 

RESULTS AND DISCUSSION 

Figure 4 shows mean streamlines for the SF, SB, AD, and computed cases of a = 5°. 
Figure 5 presents the same information for a = 23°. The streamlines are constructed 
from the appropriate ensemble-averaged velocity (u and u) field using a 4-point spatial 
interpolation, predictor-corrector scheme. As the computations predict flows that are 
completely steady, the corresponding streamlines also represent streaklines and parti- 
cle paths. The sharp turns near the surface of the flare in the 5° streamlines reflect a 
scarcity of data in that region. The double core in the 23° AD case was observed to 
disappear when a small change in vertical velocity was imposed, and is not thought to 
accurately represent the flow. The authors believe that the overall smoothness of the 
streamlines, when no special measures were taken to achieve it, indicates the data to 
be of high quality. 

A significant point to be made from these figures is that in the experiment, there is 
little, if any change in the appearance of the mean velocity field as the shock wave 
fluctuates between its upstream and downstream positions. The recirculation core 
(that point, or region, shown by a confluence of streamlines, from which the flow exits 
the symmetry plane) moves back and forth along the cylinder. Its strength seems to 
vary also, with more of the surrounding fluid affected in the shock forward case than 
in the shock back. However, the overall form of the flow apparently does not change. 

The computations, on the other hand, depict a flow that is significantly different from 
any of the measured cases for either flare angle. The predicted recirculation zone is 
larger than even that shown in the SF case, with the particles apparently spiraling 
inward much more slowly toward a core which is located well up along the flare rather 
than over the cylinder. 

In figures 6 and 7, contours of the mean streamwise velocity are shown for the sepa- 
rated regions of q = 5 C and a = 23° respectively. Only experimental results (SF, SB, 
and AD) are given. The contours show the separation “bubble” to expand and con- 
tract like a ballon (rather than translating up- and downstream, for example) as the 
shock position varies. Relative to the motion of the boundary-layer separation point 
(approximated as where the u = 0 contour is nearest the cylinder, and marked “S”), 
reattachment (u = 0 nearest the flare, marked “R”) seems to vary far less, from SF to 

SB, for a = 23° than for a = 5°. This may be another indication of the protuberance- 
like “character of the 23° flare which was discussed previously in some detail. 6 

Figures 8 and 9 show the streamwise distributions (SF, SB, AD, and computations) 
of the maximum turbulent kinetic energy ( k max ) for the two flare angles studied. It 
is noted that the experimental kinetic energies are from 2-component velocity data, 
and the assumption has been made that the third normal stress, w' rmi is equal to the 
average of the two that were measured. 

For a = 5°, it is almost uniformly true that k max for the SF case is larger than kmax 
for AD at a given location. The SB value is consistently the lowest. The reasons for 
this are presently unknown, but the digitization of high-speed shadowgraph movies, 
and the computer simulation of a moving shock w r ave are being used to search for 
clues. 



x, cm 



Figure 4 Mean streamlines, a = 5°, <f> — 0°. 

The fact that k max for the AD case exceeds the SF value at x = —3 cm, is due to 
strong influence from the shock motion. The distribution of velocities at this point ex- 
hibits two distinct states, such as Si and S? in the joint probability density diagram 
shown in fig. 10. Each state can be assigned its own probability of occurrence, Pi and 
P 2 , along with values for mean velocities and turbulent stresses (i.e., uj, V 2 Tmt , etc.). 



Reference 7 asserts that the turbulent kinetic energy of the all-data case is not merely 


a weighted (by P) sum of the energies for the two states which comprise it, but is also 
a function of the difference between the states and the product of their probabilities. 
This “organized” or “coherent” (as opposed to dissipative) contribution to the turbu- 


lent kinetic energy is expressed by PiP 2 ((^i — U 2 ) 2 + (vj — t^) 2 )* It is clearly significant 
in the immediate vicinity of the shock, where histograms such as that shown in fig. 10 


are common, but its importance appears to fall off rapidly away from the shock. 


The distribution of k max for a = 23° (fig. 9) is more complicated than for 5°. The AD 
levels remain higher than those for SF from x = —5.8 cm to i = —0.3 cm. This may 
result from larger amplitude shock oscillations, and a more persistent “coherent” con- 
tribution. However, for much of this same span, the SF levels stay at or below those 
for the SB case. The trends for the 5° flare are closely matched by those observed for 
o = 0° and 10° (though the latter two are not shown), and once again the 23° flow 
demonstrates significantly different behavior than for flows over lesser inclined flares. 
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Figure 10 Joint probability density diagram from LDV data, 
4> = 0°, a = 5°, x = 3 cm. 


Finally, the computed k max distributions bear little or no resemblance to those mea- 
sured, whether conditionally sampled or not. Even when compared to the SB case, 
for both flares, the initial rise in computed turbulent energy occurs far downstream of 
where it is actually measured. The absolute peak k levels in the computations are w'ell 
below 60%, of those measured, and they also occur at different x-locations. 

CONCLUSIONS 

The authors have studied numerous two- and three-dimensional shock-wave turbu- 
lent boundary-layer interactions at Mach 2.85. This paper has presented results for 
two fully separated flows, one moderately three dimensional, and the other highly so. 
The analysis focussed on the large scale unsteadiness of the separation shock wave 
(which has long been observed in this class of flow's) as it relates to time-averaged 
flow quantities such as mean velocities and turbulent stresses. To do so, a basic tech- 
nique was employed to conditionally sample flow-field measurements using an instan- 
taneous wall pressure level, and statistics calculated from the resulting sorted ensem- 
bles were compared with those derived from the unsampled data, and with computed 
(Navier-Stokes) results. Some basic questions, outlined in the introduction, have been 
addressed. 

Regarding the ability to remove unsteadiness effects through conditional sampling: 

The coupling of high speed shadowgraph and surface-pressure data indicated that the 
latter could be used to infer shock position. The performance herein of a relatively 
simple conditional sampling scheme which uses an instantaneous pressure level to sort 

velocity data, suggests that “snapshots” of the flow corresponding to specific shock 
positions can be obtained. This, however, does not remove all unsteadiness effects, 
since the many fluid particles whose velocities comprise those snapshots still were in- 
fluenced by an unsteady shock wave. 


On the subject of “pseudo-turbulence”: It is possible that for the a = 23° case, which 
is highly three dimensional, turbulence quantities are dominated by a “coherent” 
contribution from the shock motion throughout much of the flow. For the lesser in- 
clined flare, any such effect seems to be limited to the immediate vicinity of the shock, 
the only place where the all-data case shows higher turbulence levels than the sorted 
shock-forward (SF) case. 

Finally, on the issue of how important the shock unsteadiness is to the accuracy of the 
computations, the following points should be considered. As the shock moves forward, 
the reverse flow region seems to expand, pushing the separation point forward along 
the cylinder, the reattachment point aft along the flare, and the zero velocity line out- 
ward from the corner. The region contracts in a similar fashion as the shock retreats. 
Still, the basic character of the mean flow-field, identified by a recirculation core lo- 
cated over the cylinder, does not change with shock position. However, the computa- 
tions predict a completely steady flow, with the recirculation core above the flare, that 
does not approximate that basic character. This gives cause to wonder whether any 
improvement in computational accuracy would be realized simply by the addition of 
an unsteady shock wave mechanism. 
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